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The L3 data on Bose-Einstein correlations of equally charged pion pairs 
produced in hadronic Z decays are analyzed in terms of various parametriza- 
tions. Preliminary results are presented here. 
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i ' 1. Introduction 

QJ , In particle and nuclear physics intensity interferometry provides a direct 

experimental method for the determination of sizes, shapes and lifetimes 
of particle emitting sources (for recent reviews see |2]). In particular, 
boson interferometry provides a powerful tool for the investigation of the 
■ space-time structure of particle production processes, since Bose-Einstein 

correlations (BEC) of two identical bosons reflect both geometrical and 
dynamical properties of the particle radiating source. 

Originally, the method of Bose-Einstein correlations was invented by the 
radio astronomers R. Hanbury Brown and R. Q. Twiss (HBT), who applied 
it to determine the angular diameters of main sequence stars [31 H] . The first 
experimental evidence for Bose-Einstein correlations in high energy physics 
dates back to 1959 The results were interpreted as BEC by G. and S. 
Goldhaber, Lee and Pais (GGLP) 6 j. Angular distributions of pions could 
be described more exactly by applying Bose-Einstein statistics instead of 
a classical statistical model. This effect is frequently referred as either the 
HBT, or GGLP effect, or simply Bose-Einstein correlations. 

The L3 experiment provides a good opportunity to perform detailed 
investigation of Bose-Einstein correlations in e + e _ annihilation at a center 



* Presented by T. Novak at the 4 th Budapest Winter School on Heavy Ion Physics, 
December 1 - 3, 2004. To be published in Acta Physica Hungarica Heavy Ion Physics. 



(1) 



2 



ARXIV-RHIC04 PRINTED ON FEBRUARY 5, 2008 



of mass energy of y^i = 91 GeV. After "standard" event and track selection 
we concentrate on 2-jet events detected using the Durham algorithm [7] 
with resolution parameter y cu t = 0.006. Finally, approximately 13 million 
pion pairs are analyzed. The results presented here are preliminary. 



2. Shape of the BEC function 

The two-particle Bose-Einstein correlation function is defined as: 

C 2 (P1,P2) = - ( r— ( (1) 

where P2(pi,P2) is the two-particle invariant momentum distribution, p\{pi) 
is the single-particle invariant momentum distribution and pi stands for the 
four-momentum of particle i. Since it is difficult to create the product of 
the single-particle distribution it is replaced by a "reference sample", the 
two-particle density that would occur in the absence of BE interference. 

If long-range correlations can be neglected or corrected for, and if the 
short-range correlations are dominated by Bose-Einstein correlations, this 
two-particle BEC function is related to the Fourier-transformed source dis- 
tribution. If we assume that f(x) is the density distribution of the source 
of the pions then the correlation function is found to be 

C2(pi,P2) = l + \f(Q)\ 2 , (2) 

where Q is the invariant four-momentum difference, Q = — {pi — P2) 2 and 
f{Q) is the Fourier transform of f{x) 

f(Q) = J dxexp(iQx)/(x). (3) 



2.1. Gaussian distributed source 

The simplest assumption is that the source has a Gaussian distribu- 
tion, in which case the Fourier transformed source function is determined 



as f(Q) = ex P ( — which yields 



C 2 (Q)=7 



1 + Aexp (-R 2 Q 2 )] {1 + SQ) 



(4) 



where the parameter 7 is a constant of normalization, A is an incoherence 
factor, which measures the strength of the correlation, R is a scale parameter 
which measures the length of homogeneity and (1 + SQ) is introduced to 
parametrize possible long-range correlations. 
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Fig. 1. The Bose-Einstein correlation function C2. The full line corresponds to the 
fit by Eq0] while the dashed line stands for the long-range momentum correlations. 

A fit of this correlation function results in an unacceptably low confi- 
dence level from which one concludes that the shape of the source deviates 
from a Gaussian. The fit is particularly bad at low Q values, as shown in 
Fig. □ 

2.2. Levy distributed source 

The Central limit theorem states that, under certain conditions, the sum 
of a large number of random variables behaves as a Gaussian distribution. 
The generalization of the central limit theorem gives the classification of the 
stable distributions. The study of these stable distributions was begun by 
Paul Levy in the 1920's. For recent results see jSJ. 

According to the generalized Central Limit Theorem, if particle produc- 
tion is a result of a multifold probabilistic process such as the branching 
of gluons into gluons into gluons etc., then the class of possible limiting 
probability distributions coincides with the class of Levy distributions |S] • 
The characteristic function of symmetric stable distributions is jH] 



f(Q) = exp (iQ8 - i^L!) 



(5) 
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Here we utilize the notation of Nolan 8/,. The index of stability, a, satisfies 
the inequality < a < 2. The case a = 2 corresponds to a Gaussian source 
distribution. For more details see |S]. 

Thus the Bose-Einstein correlation function for Levy stable distributions 
has the following, relatively simple form: 



C 2 (Q)= 7 [l + Xexp(-{RQ) a )](l + 5Q) 



(6) 



After fitting the data with Eq|S]it is clear that the correlation function is far 
from Gaussian: a ~ 1.3. However, the confidence level is still unacceptably 
low. The curve of this fit is shown in Fig. |^J 
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Fig. 2. The Bose-Einstein correlation function C2. The full line corresponds to the 
fit by Eq0J while the dashed line stands for the long-range momentum correlations. 



Since there is no particle production before the onset of the collision, a 
more appropriate form of the source distribution for the time component is 
the asymmetric stable distribution. In this case, one obtains the following 
result for the correlation function: 



C 2 {Q) = 7 [1 + Acos [{R a Q) a ] exp {-{RQ) a )} (1 + 5Q), 



(7) 



where R a is an additional parameter. 

The fit of Eq|7| to the data is statistically acceptable. The data points 
are well described by the fit curve which is shown in Fig. |21 Note that for 
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Fig. 3. The Bose-Einstein correlation function C^. The full line corresponds to the 
fit by Eq0] while the dashed line stands for the long-range momentum correlations. 

Q between 0.5 GeV and 1.5 GeV the fitted curve goes below one. This is 
caused by the cosine term, which comes from the asymmetric Levy assump- 
tion, in Eq0 The fitted value of the index of stability, a is found to be 
approximately 0.8. 



The assumption that the source has a Gaussian shape is too simple. A 
good description of the Bose-Einstein correlation function is achieved using 
Levy stable distributions as the source function. 
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